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We suggest a mechanism explaining the origin of three generations of the 
Standard Model fermions from one generation in a higher-dimensional theory. 
Four-dimensional fermions appear as zero modes trapped in the core of a topo- 
logical defect with topological number three. We discuss hierarchical pattern 
of masses and mixings which arises in these models. 



1 



1 Introduction 



Recently, the interest has been renewed to particle physics models in more than 
four spacetime dimensions (see, for instance, Refs.|l[]). These models provide an 
interesting framework for solving hierarchy problems of the Standard Model of particle 
interactions. In particular, it has been pointed out |0 that by introducing extra 
dimensions it is possible to explain the hierarchical fermionic mass pattern. The 
Standard Model fermions as well as Higgs boson are represented by localized modes 
in extra dimensions [0], and their effective four- dimensional Yukawa couplings are 
determined by overlaps between Higgs and fermionic wave functions. To produce 
hierarchical structure in overlaps, hence in four-dimensional Yukawa couplings, it has 
been proposed to localize Higgs and three fermionic generations at different points in 
extra dimensions 

In this paper, we follow a completely different approach to produce fermionic spec- 
trum via wave function overlaps, explaining simultaneously the origin of three gener- 
ations of fermions with similar quantum numbers. Suppose one has single fermionic 
generation in a multi-dimensional theory. Let us consider a topological defect whose 
core corresponds to our four-dimensional world. Chiral fermionic zero modes may 
be trapped in the core due to specific interaction with the fields which build up the 
defect 1^, 1^, ^ . In some cases, the index theorem guarantees that the number of chiral 
zero modes is equal to the topological number of the defect (see for instance Refs.Q]). 
We will use this property to obtain three fermionic generations localized on a defect 
with topological number three while having only one generation in the bulk. If the 
Higgs scalar couples to the defect, it can also be trapped in the core Hierarchy 
between masses of three fermionic modes arises due to their different profiles in extra 
dimensions. 

To be specific, we will work with the simplest topological defect which can be 
characterized by an integer topological number, the global vortex (see Sec. 2). We 
first recall the mechanism of localization of fermions (Sec. 3) and scalars (Sec. 4), 
then explain the origin of different Yukawa couplings for similar fermions of three 
generations (Sec. 5). To explain mixing angles between fermions of different genera- 
tions, we then introduce a complication in the model (Sec. 6). We present a general 
discussion of the presented class of models in Sec. 7. Notations and technical details 
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are outlined in Appendices. 



2 A global vortex with topological number k. 

Consider a theory of a complex scalar field $ in six dimensions, 

u = \d^^f-^{m'-v')\ (1) 

(See Appendix ^ for notations). The global U{l)g symmetry $ — > ^e*"" is broken 
spontaneously by the vacuum expectation value |$| = v. 

Let us consider field configurations which do not depend on x^. We introduce 
polar coordinates (r, 9) in (0:4, 0:5) plane. There are solutions to the classical field 
equations which have the form 

$ = ye*'=^F(r), k = ±l,±2,.... (2) 

Their topological numbers are defined by the winding numbers k; the function F{r) 
satisfies the following boundary conditions: 

F(r) ^ 1, r ^ 00; 

F{r) 0,r 0, 
and the ordinary differential equation, 

F" + -F' -^F- Xv^F (F' - 1) = 

(hereafter prime denotes the derivative with respect to r; we assume > 0). Ana- 
lytical solution to this equation is unknown. However, from this equation, it follows 
that 

F{r) = 0{r^), r ^ 0; 
F{r) = l-^ + 0{r-'), r^oo. 
So, the field configuration describes a "3-brane" in six-dimensional spacetime, 
located at X4 = X5 = 0, of radial size of order {^\/Xv^ 

Energy of this configuration diverges logarithmically at large distances, r 00. 
This can be cured either by gauging U{l)g (the case of local, or Abrikosov-Nielsen- 
Olesen vortex) or by introducing some cutoff at large distances (for example, putting 
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an "anti-vortex" with winding number {—k) far away from the vortex). In any case, 
this would improve the behavior of energy density at large distances. Here, we are 
interested in physics inside the core which remains intact after these modifications. 

3 Fermionic modes. 

Consider a six-dimensional fermionic field Q which has an axial charge 1/2 under 



We will assume that g is real: its complex phase can be rotated away by an axial 
transformation of the fermion. 

Let us perform Fourier transform with respect to four- dimensional coordinates x'^, 



U{1) 



and interacts with $ via axial Yukawa coupling. 





The Dirac equation following from the Lagrangian Eq.(^) takes the form of Schrodinger 
type equation for \l/ in the vortex background (^, 



(5) 



where C' = T^r, and 




The operators Cj and D anticommute. 



QD + DCi = 0. 



This means in particular that one can search for a solution of Eq. as an expansion 
in the eigenvectors of the operator D, 



m 



m 



(6) 
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There may exist a set of discrete eigenvalues m with separation of order gv, and 
continuous spectrum starting from m > gv. For given A;^, the eigenvalues m satisfy 
^0 ~ + This means that to excite a mode with non-zero m, energy of order 
gv is required. In what follows we assume that the energy scales probed by a four- 
dimensional observer are much smaller than gv, and thus even the first non-zero level 
is not excited. So, we are interested only in zero modes of D: 

= . (7) 

As is shown in Appendix 0, the solution for Q{x^, X4, xs) has the form 



p 



( \ 

/p(r)e^f^Ck 
/fc-p-i(r)e-*('=-P-i)Vk 

V / 



(8) 



where Ck is a normalized four-dimensional left-handed spinor, are arbitrary complex 
functions of a three-dimensional momentum k = {/ci}, and ]i{r) is the normalized 
solution to the following differential equation, 



At r — > 0, /z(r) ~ r', while at r ^ 00, fi{r) ~ e"^''^. 

In other words, there exist k zero fermion modes with different angular and radial 
wave functions. One bulk fermion corresponds to k four-dimensional fermion species. 

Let us emphasize that at given p, there is only one degree of freedom for particle 
(fco > 0) and one degree of freedom for antiparticle (/cq < 0) that corresponds to one 
spin state of a left-handed spinor. Since /p(r) fall off exponentially at large r, the 
particles described by are localized in the core of the vortex. 

It is instructive to study how this eight-component localized spinor can be con- 
sidered from the four-dimensional point of view. For this purpose, let us introduce a 
8x8 representation of four- dimensional 7-matrices: 7^ = F'^, 



7' = zFoFiFaFg 
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These matrices together with the unit matrix, their commutators a^'^ = |[7'^7'^], 
and 7^7^ form an algebra A (that is any product of these matrices is their hnear 
combination). The algebra A is isomorphic to algebra ^ of 7^^, and one can construct 
an operator which relates 7^^ and 7^^: 



t/+7t/ = 7 



(10) 



for any 7 G 7 G ^. One can check that 



U 



1 / 7o(l - 75) 
1 + 75 



/O l\ 



1 
\0 0/ 



:ii) 



is one of the solutions to Eg. (pUj) . To find four-component spinor uj{x) which corre- 
sponds to fl{x), one should act by on ^2: 



re 



ip0 



Ck 





fc-i 

p=0 

We again obtain the same result: u describes k left-handed massless spinors ^Jpi^x^) 
localized at the vortex. 

To localize k right-handed massless spinors at the vortex, one should consider 
six-dimensional spinor S which has an axial charge —1/2 under U{l)g: with obvious 
modifications, the above analysis goes through for S as well. 



4 Higgs field. 

Let us introduce a complex scalar field H which interacts with the vortex field 



H 



\dAH\^--{\H? 



K 

2 



(12) 



Note that at this point H can have arbitrary charge under U{l)g. Later on, we will 
set this charge to zero. 



6 



The system of the two scalar fields, $ and H, described by the Lagrangian C = 
jC<i, + Ch, Eqs.(|I]) and (|12|), admits several nontrivial classical solutions. One of them 
corresponds to $ given by Eq.(^ and H = 0. It can be shown 0, however, that 
this solution is unstable in a certain region of parameter space, in particular, for 
KyU^ < h'^v'^. We will consider this case in what follows; the lowest energy solution in 
the topological sector where $ has the form (|^) is 

where radial functions Fc{r), Hc{r) satisfy the following set of nonlinear differential 
equations, 

K' + -rK - S^c - h-'FMl - Xv'F, (F2 - 1) = 0, 
H'^ + \H'^ - h^H.v^F^ - kH,{H^ - fi^) = 0. 

The boundary conditions are 



(14) 



F,(0) = 0, F,(oo) = l; 

where vh is the vacuum expectation value of the Higgs field from four-dimensional 
point of view. The leading behavior of Hc{r) near the origin is 

H{r) = vn(l- <ll^A . (15) 



So, to satisfy to the boundary condition at r — > oo, it is required that vh < /i- The 
leading behavior of Fc{r) is the same as of F{r), Eq.(|^). In the background of the 
modified vortex solution, Eq.(|13|), one still recovers all results of Sec. |], replacing 
F{r) by F,(r). 



5 Four-dimensional fermion masses. 

Let us turn to the fermion-Higgs Yukawa couplings. To describe a "prototype" gen- 
eration, we need two six-dimensional spinors, Q and U, with opposite charges under 
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The interaction with the vortex with k = 3, 

^,$gi^g + g^<i>*u^—^u + h.c. 

results in the existence of three left- and right-handed zero modes of Q and U, 
respectively: 

/ \ 



Pq=0 



Ckg2-p,(r)e~^(2-P5) 



v 







(16) 



U 



E 

Pu=0 



(27r)3/2 



j rf'A;5(A;2)e-*'='*"^4^)(k) 








Here, Ck and dk are the left-handed and right-handed four-dimensional spinors, respec- 
tively; arbitrary complex functions al?](k) for each pq = 0, 1,2 describe two degrees 
of freedom (one for the particle and one for the antiparticle) of a massless four- 
dimensional left-handed fermion Qp^ with gauge and global quantum numbers of Q; 
a"p] (k) describe a massless four-dimensional right-handed fermion Up^ with quantum 
numbers of U\ Pq = 0, 1, 2 corresponds to three generations of fermions. The wave 
function profiles in extra dimensions are defined by radial functions qpg{r), Up^{r), 
which are the solutions of Eqs.(^ with k = 3, f = q or u, g = gq or gu, and I = Pq or 
Pu, respectively. 

Fermion masses originate from the following term in the six-dimensional action: 



t/ + h.c. 



{Yu is the six-dimensional Yukawa coupling constant, and Hi = eijH*). This inter- 
action is U{l)g invariant only if H is neutral under U{l)g. To obtain the effective 
four-dimensional mass matrix, one has to perform the integration over extra dimen- 
sions in the action integral. This results in the Dirac mass terms. 



m 



QpqUp 



PqPu^Pq'-^Pu 1 
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with 

= YuJrdrdeH,{r)qp^{r)up^{r)e'^P''-P-^'. (17) 

Integration over 6 in Eq.(|l^ leads to a selection rule, rrip^p^ ~ ^PqPu- This means 
that in this way it is possible to generate only diagonal mass terms, but not inter- 
generation mixings. To obtain mixings between fermions of different generations, it 
is necessary to have a non-trivial 6 dependence in the Higgs mode. We will return to 
this point in Sec. 6. 

To see that the hierarchy of diagonal mass values can be generated in this way, one 
can make use of the following, very rough, estimation. Characteristic distance scales 
for the Higgs and fermionic modes are {^/K,f^)~^ ~ {hv)~^ and {gv)~^, respectively. 
Let us take h ^ g, so that the Higgs mode is narrow in comparison to fermionic 
modes. Then one can substitute fermionic radial functions in Eq. ( |T7|) by their leading 
behavior at r — 0, 

Qp ~ {gqvry, Up ~ {guVrY 
and to write He = Hc{hvr). In this approximation, one gets 

(18) 

up to some slowly varying function on p,g, and h. Thus, at 5 ~ 0.1, mass hierarchy 
among three generations is naturally reproduced, if we associate the first generation 
with p = 2, the second one with p = 1 and the third one with p = 0. To obtain exact 
predictions, one has to solve the equations (^, (p^ and to evaluate the radial integral 
in ( p!7D numerically. A similar procedure yields masses of down-quarks and charged 
leptons, starting from interactions Eqs. (0), (0) (see below). Since dependence 
on g of wave function overlaps is highly nonlinear, the model does not predict any 
simple analytic relation between masses of fermions of different kinds (for example, 
mt/rric ^ mb/rris, etc.). 

The parameters of the model with all Standard Model fermions included are: 
an overall mass scale, say, y^y^/i; two ratios of six- dimensional Yukawa couplings, 
Yd/Yu and Yi/Y^, and five ratios gq,u,d,i,e/h. As follows from the discussion above (see 
Eq. (0)), to the leading approximation, masses depend only on three combinations 
of the latter five ratios, namely, 5„ = gggu/h"^, = gggd/h"^, and 6i = gige/h^- Thus, 
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diagonal masses of nine charged fermions of the Standard Model arise from eight 
independent parameters, with significant dependence from six of them. 

6 Mixing between generations. 

To obtain off-diagonal mass matrix elements which mix fermions of different genera- 
tions, one has to relax the selection rule ^p,p„ in Eq.(|l^) by introducing non-trivial 9 
dependence in the Higgs mode. To do this, one has to complicate the model, because 
the interaction (^2|) is phase-independent, so the classical solution He depends only 
on r. In what follows, we will need a somewhat more complicated topological defect, a 
global vortex made of two scalar fields with different winding numbers which appears 
in a model with 

C = jC^ + Cx, 
where is presented in Eq.(|I]), and 

Complex six- dimensional scalar fields $ and X have charges 3 and 1 under global 
U{l)g symmetry, respectively^. Let v and vi be vacuum expectation values of |$| and 
|X|, respectively; 

V = V + 0(a), 
vi = vi + 0{a). 

At a < Ai, the model admits stable global vortices of the form 

(19) 

X = fie x(^); 

with F{r) different than F{r) in Eq.(H), but still having the same leading behavior 
at r — > and r oo. Note that in this model, the solution, Eq. (|19|), is the simplest 
nontrivial vortex. 

Let us couple six-dimensional fermions corresponding to five fermionic fields of 

one generation of the Standard Model (namely, left-handed leptons L and quarks Q 

*We assume that a is real; this can always be reached by rotating away its complex phase via 
redefinition of the origin of polar angle 6. 
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and right-handed leptons E, up quarks U and down quarks D) to $ in a way similar 
to Eq.(|), 

gu^*U^U + gd^*D^D + ge^*E^E + \i.c., 

where we have written exphcitly the electroweak SU{2) indices i, and to the Higgs 
field H, 

Vuf = Y^H,Q,^—^U + YaH,Q,^-^D + YiH,U^-^E + h.c, (21) 

where Hi = €ijHj. Three chiral generations are localized in four dimensions as zero 
modes of the form ([T6| ) (with F{r) replaced by F{r) in corresponding equations). 

On the other hand, let the Higgs field H couple to the field X in a way similar to 
Eq-dH), 

V^H = ^im'- f^'Y + hW\X\' + AV, (22) 

where AV is a small perturbation which breaks U{l)g global symmetry explicitly, for 
instance, 

AV^ = \Hi\'^{eX + e*X*). 

At e = 0, we recover Eq.(|I^ with Fc(r) replaced by x(^); its classical solution, call 
it Hglr), contributes to the diagonal masses through the integral Eq.(|l^. If e 7^ 0, 
the solution to the nonlinear partial differential equation for H{r, 9) which follows 
from ( P^ ) in the background (|I^) has a complicated dependence. For our purposes, 
however, it is sufficient to perform perturbative expansion, 

H{r, 6) = Ho{r) + eH,{r, 9) + e^H^ir, 9) + ... + h.c. (23) 

The first term is Ho{r) defined above. In the first order in e, a source term of the 
form 

appears in the equation for Hi{r, 9) which results in 

Hi{r,9) =e'^hi{r). 

In a similar way, 

H2{r,9)=e^''h2{r). 
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The integral (|T^) reads nowQ 

m^f, =YuJrdrde H*{r, e)g3_a(r)u3-fe(r)e*(''-'')^ = 

= Y^,{2tt) Jrdrq3_a{r)u3_b{r) {HQ{r)6ab+ ^^4) 

e/ii(r)5a-i,fe + ^*hiir)6a,b-i+ 
e^h2{r)5a-2,b + t*^h2{r)5a,b-2 + ••.)• 

Non-diagonal mass matrix elements arise from overlaps of /ii,2(^) with fermionic wave 
functions and are suppressed by powers of e with respect to diagonal ones. 

In the form discussed above, the model does not admit CP violation since the 
overall phase of is irrelevant and phase of e can be included in two-component 
spinors Ck and (see Eq.(0)). With more complicated A\^, for example, 

= eX\Hi\^ + eiX'^\Hi\^ + \i.c., (25) 

the relative phase of ei and e is a free CP violating parameter of the model (while all 
other results remain intact). 

The relevant real parameters are now gq,u,d,i,ev/ (hvi), |e|, |ei|, arg(ei — e) and 
again Y^^/HyL, and Y^i^ilYu. These eleven parameters of the six- dimensional theory 
thus generate thirteen parameters of the fermionic sector of the Standard Model (nine 
diagonal masses, three mixing angles, and one CP violating phase). As before, results 
depend significantly not on all five ratios gv/{hvi), but only on three combinations 

ggguV^ihW), gggdV^ihW), gWev'^/ih^vi^)- 

To demonstrate that it is possible to obtain a realistic hierarchical pattern of the 
fourteen Standard Model masses and mixings, let us estimate them in more detail. 
First of all, we find H{r, 6) (see Eg. (^31) ) in the model with AV given by Eq.(^). We 
assume that e and ei are of the same order and find H{r, 6) in the first order in e, ei: 

H{r, 6) = Hoir) + h{r){ee'' + h.c.) + h2ir){e,e^'' + h.c). (26) 

Close to the origin 

Ho{r) = Vh] h{r) = r^; /i2(r) = — r , (27) 



^Hereafter, the indices a,b = 1,2,3 are used instead oipq,Pu,--- = 2,1,0 to enumerate the 
fermionic generations; a = 1 corresponds to the first generation, etc. 
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(recall that X has winding number one and thus X'{0) ^ 0). On the other hand, the 
width of hi^2 is of order of width of Hq that is of order (hv)^^. Substituting Eq.(p7D 
into Eq.(^) and estimating all integrals as it has been done it in Sec. ^ we obtain 
the following mass matrix of up-quarks 



ab 



Yu 



3 \ 



e6i 5u e*5i 
\ ei^n 1 j 



and a similar matrix with the replacement 5u —>■ S^, —>■ Y^i, for down quarks. 
Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix is defined by 



U 



CKM 



S},Sd 



where Su and Sd transform the matrices mum\^ and mdm\^ respectively, to the diag- 
onal form: 

S^mm^ S = diag(. . .). 

With the definitions 



^21^11 + ^12^^22 



m 



u2 
22 



m 



u2 
11 



~ e 



^23^33 + 



22 



7/9 1/9 

— 



7/9 7/9 

777,33 ~ "^11 



J 



and the same for down quarks we find that in the leading order, 

^ 1 (Ad-A^) (C,-a) ^ 

f/™^ (Au-Ad) 1 (i?,-i?Je-^ 

^ {Cu - Cd) {Bu - B^y^ 1 

I 1 mi-^i) \eim-^i) \ 

where (/? = argei — 2arge. We see that there is a hierarchy in CKM matrix in our 
model which coincides with the hierarchy of CKM in Standard Model: the magnitude 
of the matrix elements decreases away from the diagonal, so that Vcb ~ Kjs ~ 6(5^, 
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7 Conclusions. 



The mechanism discussed here provides a consistent picture which can explain the ori- 
gin of three generations of fermions with identical gauge and global quantum numbers 
but hierarchical mass matrices without fine tuning of parameters. Chiral fermionic 
zero modes are localized on a topological defect with topological number three which 
explains the origin of three generations; different masses appear due to different pro- 
files of three fermionic zero modes in extra dimensions. 

To explain mixings between generations, one has to introduce small but explicit 
violation of the global U{l)g symmetry of the theory. Off-diagonal mass matrix 
elements are suppressed by powers of the small parameter e characterizing this viola- 
tion. Though topological arguments are no longer valid in the presence of this U{l)g 
violation, we expect our results to be stable with respect to e. 

The model with a global vortex can be embedded either in a theory of large 
compact extra dimensions (in which case the problem of stability of the vortex in a 
compact space has to be addressed) or in a model with localized gravity [0]. The 
mechanism suggested here works as well in models with other topological defects 
with topological number three, for example, with a "hedgehog" in seven space-time 
dimensions. 

The models which exploit the mechanism discussed above are in principle fairly 
predictive. For instance, in our particular toy model with a global vortex, eleven 
independent parameters of the six-dimensional theory generate thirteen masses and 
mixings of the Standard Model fermions, and the hierarchical structure of mass ma- 
trices is reproduced. 
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A Notations. 

Six-dimensional coordinates xa are labeled by capital Latin indices A,B — 0, . . . , 5. 
Four-dimensional coordinates are labeled by Greek indices n,!/ — 0, . . . , 3; for 
spatial coordinates we use lower case Latin indices i,j = 1,2,3. The Minkowski 
metric is qab = diag(-F, -,...,-). 

Dirac fermions in six dimensions are described by eight-component spinors; we 
work with the following representation of six-dimensional 8x8 Dirac matrices F"^: 




where = T.^ = 7V; = = 7V; = -S^ = i7°7^ = -E^ = 7°, and 
7^*, 7^ are the usual four-dimensional Dirac matrices in the chiral representation: 

7=1 I, ^ - qI' ^ =^0717273 = 

(7* are the Pauli matrices. 

We also introduce Fy which is an analog of four-dimensional matrix 75 



r7 = ro...r5 = 
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B Fermionic zero modes and their asymptotics. 

The solution to Eq.(|^ has the following structure 



( /(i)(r)c(i)e^(P+i)^ \ 
/(2)(r)c(2)e^^'^ 
/(3)(r)c(3)e-(^-i-P)^ 
V /(4)(r)c(4)e-*('^-f)^ / 



(28) 



In the last equation, p is an integer number; C(a) are two component columns (which, 
as will be shown below, correspond to four- dimensional chiral spinors), and /(q) satisfy 
the following set of differential equations |, 



I /(4) + ^^/(4) - gvFf^,^ = 0; 
' /(2) - zh) + 9vFU,) = 0, 



(29) 



(30) 



To investigate the behavior of /(„) and to find their asymptotics, it is convenient 
to introduce the new set of functions /(„): 

/(I) = r-(^+'V~(i) , h)=r'~h), fi3) = r^'-'-'^fi3), /(4) = r-('=-^)/(4) • (31) 

The functions f\a) satisfy the following set of differential equations obtained from 

■l-These equations cannot be solved analytically, except for a particular case oi p — {k — l)/2, 
when a normalized solution has the form 



/(2),(3) = const • r 2 exp - / gvF{x)dx , /(i),(4) = . 
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Eqs.(Pi, (Pi: 



/(4) - 

/(I) = 

/(2) - 

/(3) = 



F' {k-2p- 1] 
T ^ r 



/('4)-(7VFV(4) = 0, 



(32) 



(4) 



F' {k-2p-l] 
f ^ r 



f[2) - 9'v'F''fi2) = 0, 



(33) 



f. 



(2) 



Note that /(i) and /(3) satisfy the same differential equations as /(4) and /(2) 
respectively, with p replaced by A; — p — 1. 

/^From Eqs. (|32|) it follows that, if at some point tq a solution /(4) and its derivative 
are positive (/(4)(ro) > 0, /(4)(ro) > 0), then /(4) increases with r at any point r > tq. 
To see this, let us first note that equation on /(„) is nothing but a Schrodinger type 
equation for the lowest energy level. This means that /(q) cannot have nodes at 
< r < oo, so we can assume /(4) > 0. Let us assume also that (fc — l)/2 > p, and 
thus k>p. Since F' / F > and /(4)(ro) > 0, it follows from Eq.(|3|) that /(4)(ro) > 0. 
It means that /^'^-j increases and so /(4)(t') > at any point r > vq. 

Let us now study the behavior of /(a) at r — > 0. In this limit, F and from 
Eqs. (|52|) and it follows that 

/(4) ~ r° or /(4) ~ r^^'^-f) ; /(i) - r° , 

/(2) ~ r° or /(2) ~ r'^^-^') ; /(g) ^ r° . 

Since k > p, the solutions which have the behavior /(2),(4) — r^^^^^^^ increase every- 
where. 

Now let us consider asymptotics of /(„) at r — >• oo, when F ~ 1 and first equation 
in (521) reads 



/(4) /(4) -9 V f(4)=0. 

This equation has two linearly independent solutions, 

/i) = r^'-'^^^'K^.^^Mavr) ; /£ = r(^-2^)/2/(fc_2p)/2(^7^^r) , 
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where Iu{f) and K^{r) are modified Bessel functions of tlie first and tlie tliird order. 
Tlie second solution /^Y) grows exponentially at infinity. On the other hand, as we 



have shown above, if a solution to equations (P^), (33) increases at some r, then it 
increases everywhere. This means that the solution which behaves as /(4) ~ r'^i^~^'> 
at r — has the asymptotic behavior /^Y) at r — > oo. Thus, the corresponding 
function is not normalizable. Therefore, all functions /(a) which can correspond to 
normalizable solutions have the same asymptotics: they tend to constant at r = and 
exponentially fall off at infinity. This results in the following behavior of functions 
/(a): they exponentially fall off at infinity. 



/(a) 



oc e 



-gvr 



OO 



and close to the origin, 

/(,) ^ r-(^'+^) , /(4) 



-ik-p) 



f{2) 



fi3) 



. 



(34) 



(35) 



Since k — p>0, /(4) is not normalizable: the corresponding integral diverges at r = 0. 
So, one is forced to conclude that /(i) = /(4) = 0. Moreover, to allow /(2) to be 
normalizable one should require that p > 0. 

The case 2p > k — 1 can be studied in the same way: instead of functions /(4) and 
/(2), one should consider functions /(i) and /(3). As a result one finds that k — 1 > p. 
There are k linearly independent normalizable solutions which can be labeled by index 
p = 0,...,k-l: 

( \ 



/p(r)e^PVp 



v 



/A.-i-p(r)e-('="*'-i)^c,_p_i 


where /p are the normalized solutions of the equations (|30D which behave as at 
zero and fall exponentially at infinity. 
Substituting spinor \l/p into equation 

one finds that Cp and Cfc_p satisfy the following equation 

(fco + kiai)c.p = {ko + kiai)ck-p = , 
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which defines a left-handed four-dimensional spinor for ko > 0. This equation has 
one solution (iff = kf ) which we will denote as Ck- To conclude, we have precisely 
k solutions which describe k massless left-handed four-dimensional fermions in full 
agreement with the index theorem, so the solution to Eq.(^, which corresponds to 
the zero mode, has the form (P). 
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